##### SCRIPT QUE GERA A TABELA 7 #####

# 1. CARREGA PACOTES ------------------------------------------------------

library(tidyverse)
library(janitor)
library(naniar)
library(mice)
library(xlsx)

# 2. CARREGA DADOS ------------------------------------------------------

#' Nesta etapa, foi implementada a leitura dos dados combinados de Enade-RAIS (filtrado obs com salario)

# 3. IMPLEMENTA MODELOS DE REGRESSAO OLS --------------------------------

# 3.1 Modelo geral (sem separar por areas de estudo) ----

model_sal_sem_hab <- with(dados_graduados_empregados, lm(log(salario_hora_total) ~ 
                                                               tp_sexo + 
                                                               raca_cor + 
                                                               idade_faixa + 
                                                               edu_fam))

model_sal_2 <- with(dados_graduados_empregados, lm(log(salario_hora_total) ~ 
                                                         tp_sexo + 
                                                         raca_cor + 
                                                         idade_faixa + 
                                                         edu_fam + 
                                                         # nt_ce_quartil + 
                                                         nt_fg_quartil +
                                                         trab_grad + 
                                                         setor_ies_bin + 
                                                         regiao_curso))

model_sal_3 <- with(dados_graduados_empregados, lm(log(salario_hora) ~ 
                                                         tp_sexo + 
                                                         raca_cor + 
                                                         idade_faixa + 
                                                         edu_fam + 
                                                         nt_ce_quartil_area + 
                                                         nt_fg_quartil +
                                                         trab_grad + 
                                                         igc_faixa_4c +
                                                         setor_ies_bin + 
                                                         regiao_curso))

# 4. GERA TABELAS DE RESULTADOS -------------------------------------------

combine_sal_2 <- pool(model_sal_2)
summary(combine_sal_2)
combine_sal_2_fit <- combine_sal_2$glanced

combine_sal_3 <- pool(model_sal_3)
summary(combine_sal_3)
combine_sal_3_fit <- combine_sal_3$glanced

combine_sal_3_out <- data.frame(summary(combine_sal_3)) %>% 
   mutate(coef_exp = round(exp(estimate), 3))

# 5. EXPORTA TABELA -------------------------------------------------------

xlsx::write.xlsx(combine_sal_3_out, 
                 file = "output/Tabela-7.xlsx", 
                 sheetName = "raw", 
                 row.names = F, 
                 append = T)
